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We study the effect of long-range hoppings on T c for the 2D Hubbard model with and without Hol- 
stein phonons using parameters evaluated from band structure calculations for cuprates. Employing 
the dynamical cluster approximation (DCA) with a quantum Monte Carlo (QMC) cluster solver for 
a 4-site cluster, we observe that without phonons, the long-range hoppings t! and t", generally sup- 
press T c . We argue that this trend remains valid for larger clusters. In the presence of the Holstein 
phonons, a finite t enhances T c in the under-doped region for the hole-doped system, consistent with 
LDA calculations and experiment. This is interpreted through the suppression of antiferromagnetic 
correlations and the interplay between polaronic effects and the antiferromagnetism. 
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Introduction- While most theoretical studies of the 
cuprates are in the framework of the simplest version 
of the two-dimensional Hubbard model with only near- 
est neighbor hopping, both band structure calculations 
and experimental data suggest a richer set of param- 
eters for this model iSiMiSiSiZiS. Angle-resolved pho- 
toemission spectroscopy (ARPES) plays an important 
role in this regard, suggesting different topologies of the 
Fermi surface for different high-T c superconductors 3,9 
which can be reproduced by choosing finite long-range 
hoppings (t', t", . ■ The inclusion of the next 

nearest neighbor hopping, t', in the Hubbard model is 
also necessary to capture the electron-hole asymmetry 
— iiii 2 . Furthermore, t' is an important parameter in 
determining the charge orderings and their textures in 
cuprates 



13,14,15,16 



The effect of t' on T c has been studied by different 
groups — For example, E. Pavarini et.al£ notice a 
correlation between the experimental maximum super- 
conducting temperature (T™ ax ) and the value of t' eval- 
uated from the band structure calculations in different 
cuprates. However, the mechanism which may govern 
this relationship in cuprates is not well-understood. 

Theoretical investigations employing simple models 
such as single-band Hubbard and t-J model o 12 i 18 ' 19 i 20 ' 21 , 
do not show strong evidence of a direct relationship be- 
tween T c in different doping regions and the magnitude 
of the long-range hoppings. Variety of techniques have 
been used to study the effect of t' and t" (third near- 
est neighbor hopping) on superconducting properties of 
these models. For hole-doped systems, by employing 
finite size calculations and slave-boson mean field the- 
ory, C. T. Shih et. al*£ find a strong enhancement of 
the superconducting correlations due to t' and t" in the 
intermediate- and over-doped regions and a slight sup- 
pression in the under-doped region. For electron-doped 
systems, density matrix renormalization group calcula- 
tions^ have shown that t' leads to the enhancement of d- 
wave pairing correlations. Unlike finite-size calculations 
where the transition temperature cannot be directly cal- 
culated and the superconducting properties are estimated 



approximation for the thermodynamic limit and allows 
calculation of T c and other superconducting properties 
for all doping regions. The DCA has been successful to 
derive the phase diagram of the Hubbard model, show- 
ing the antiferromagnetic (AF), pseudogap and d-wave 
superconductivity phases which are in a very good qual- 
itative agreement with experiments^^. 

In this work, we investigate the effect of t' and t" on the 
superconducting properties of the 2D single-band Hub- 
bard model with and without phonons. We find that 
without phonons, T c is generally suppressed by t' and t" . 
However, with Holstein phonons, T c increases with t' in 
the under-doped region for hole-doped systems. In other 
doping regions, phonons reduce the suppression of T c due 
to t'. The interplay between electron-phonon (EP) cou- 
pling and t' plays an essential role in the dependence of 
T c on t' and is consistent with experimental data show- 
ing evidence of strong EP interaction in the cuprate o 27 ' 28 . 
Previously, we established a synergistic relationship be- 
tween short ranged AF order and EP coupling in the 
doped Hubbard model 29 , i.e. we found that AF corre- 
lations enhance the polaronic effects (dressing of elec- 
trons by phonons) and at the same time, the EP cou- 
pling enhances the AF correlations. We also established 
that local phonons which couple to the electronic den- 
sity, strongly suppress T c due to the renormalization of 
the single-particle propagator. Here, we show that t' can 
strongly affect this synergism and thus the suppression 
of T c . A finite t' in the hole-doped systems, suppresses 
AF correlations and hence reduces the polaronic effects 
and enhances T c . 

Model- We consider a 2D Hubbard-Holstein model 
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from cluster pairing correlations, the DCA 
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is an 



where ty is the hopping matrix, c ia .(ci a ) is the creation 
(annihilation) operator for electrons on site i with spin 
cr, and U is the on-site Coulomb repulsion which is taken 
to be equal to the bandwidth (8t). We vary the filling, 
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(n), from values less than one to values larger than one 
to cover the hole-doped to the electron-doped regions re- 
spectively, ujo is the frequency of phonons and {m, Pi} are 
canonical conjugate coordinates for the phonon on site i. 
The EP coupling is on-site and proportional to the den- 
sity of electrons with the coupling strength g. We define 
the dimensionless EP coupling for Holstein phonons as 



Hole-doped 



Electron-doped 



A = g 2 /(Moj 2 Q 8t) 



(2) 



which is the ratio of the single-electron lattice deforma- 
tion energy and half of the electronic bandwidth 30 . 

Formalism- We employ the dynamical cluster approx- 
imation with a quantum Monte Carlo algorithm as the 
cluster solver. The DCA approximates the self-energy of 
the system by mapping it into a cluster of size N c em- 
bedded in a self consistent host. All of the correlations 
inside the cluster are treated non-perturbatively while 
a mean-field (MF) approximation is used to deal with 
longer range correlations. Therefore, the solution would 
be exact in the limit of N c — > oo. The Monte Carlo 
simulation performs the sum over both the discrete field 
used to decouple the Hubbard repulsion 3 -^, as well as the 
phonon field, u. Details about an efficient Monte Carlo 
simulation of systems with low energy phonons will be 
given elsewhere. 

The sign problem in QMC limits our calculations to 
relatively small clusters. Most of the calculations are 
done for a 2 x 2 cluster, the smallest cluster which allows 
wave pairing. Note that the length scale associated 
with t" is not represented for this cluster. Thus, the 
effect of t" will be similar to that in the dynamical mean 
field approximatio n 32 ! 33 ; i.e. only through changes in the 
cluster densities of states. 

In order to investigate the effect of t' on the d-wave su- 
perconductivity, we calculate the eigenvalues of the pair- 
ing matrix Txo, where xo is the bare bubble and T is the 
particle-particle irreducible vertex function calculated in 
the QMC process. At T c , the leading eigenvalue (in this 
case, the one with d-wave symmetry) goes to unity and 
causes a singularity in the two-particle pairing Green's 
function x = Xo + Xq^X = Xo/(l - ^Xo)- The value of 
the d-wave pairing interaction can be measured by cal- 
culating the d-wave projected verte x 34 i 37 
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<9(K) 2 )k 
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for the lowest Matsubara frequency and g(K) — 
cos(K. x ) — cos (K y ) where K is the momentum at the 
center of each of the N c cells which tile the Brillouin 
Zone in the DCA. To capture the effect of the dressed 
electronic propagator on the d-wave eigenvalue, we also 
calculate the d-wave projected bare bubble as 
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Results- The long-range hoppings, t' and t" can af- 
fect the superconducting phase diagram through both 
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FIG. 1: (color online) The superconducting phase diagram 
of the hole-doped (left panels) and the electron-doped (right 
panels) system for different values of t'. Upper (lower) panels 
correspond to t" — (t" = 0.2t). The open circles in c) and 
d) correspond to t' — t" — 0, and are plotted for comparison. 



the band structure and the interaction vertex. In the 
electron-doped systems, t' favors hopping in the same 
sub-lattice and enhances the AF correlations at finite 
doping while in the hole-doped systems, t' suppresses 
the AF correlationsiSi 3 ^. Presumably, t" which also in- 
troduces hopping in the same sub-lattice, would affect 
AF correlations as well. However, previous calculations 
indicate that there is a close relationship between AF 
and superconductivity in cuprate o 36 i 37 ' 38 . Therefore, t' 
and t" can influence pairing by affecting the AF correla- 
tions. They also change the band structure which alters 
the density of states at the Fermi energy and thus can 
influence T c . 

We find that t' and t" generally suppress T c in most 
doping regions, apart from a slight increase in T c at small 
dopings. First, we consider a finite t' and t" = 0. The 
superconducting phase diagrams for three different values 
of t' are shown in Fig. [T] (a) and (b). When t' = -Q.3t, 
T c is slightly smaller in comparison with the case of t' = 
from 10% to 25% hole-doping and at large electron- 
doping (>20%). With a larger t' = -OAt for the hole- 
doped system, we find that T c is strongly suppressed in 
intermediate- and over-doped regions while in the under- 
doped region, the overall effect of t' on T c is negligible. 
The effect of t" on the superconducting phase diagram is 
shown in Fig.[T](c) and (d). We see a stronger suppression 
of T c in the over-doped region and a slight increase (~ 
10%) in T c in the under-doped region when both t' and 
t" are finite. Moreover, a non-zero t" with t' = 0, does 
not have a considerable effect on T^. 

We find that the band renormalization effects due to t' 
are mostly responsible for changes in T c and the effect of 
t' on the interaction vertex is less significant. In order to 
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FIG. 2: (color online) d— wave pairing interaction, Vd (Eq.[3jl, 
and the d-wave projected bare bubble, Pdo (Eq. 2}, versus 
temperature in the under-doped and intermediate-doped re- 
gions. The empty (filled) symbols correspond to the Hubbard 
model without (with) Holstein phonons. Circles (squares) 
correspond to t' = 0.0 (i' = -0.3t). 
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FIG. 3: (color online) Inverse of the d-wave pairing suscep- 
tibility for a) 4-site and b) 16-site cluster at 15% doping for 
different next-nearest hoppings. 



illustrate this, we plot the i'-dependence of the d-wave 
bare bubble, Pdo (Eq. [4]), and the d-wave pairing inter- 
action, V d (Eq.GU, at 5% and 15% hole-doping in Fig. [3] 
(empty symbols), t' strongly suppresses Pdo and slightly 
increases Vd in both doping regions. The former effect 
is responsible for the decrease in T c . The suppression 
in Pdo is a result of the band renormalization effects of 
t 1 which decrease the density of states at the anti-nodal 
points^^. 

The effect of t 1 on the d-wave pairing shows a similar 
trend when larger clusters are considered. The inverse 
of the d-wave pairing susceptibility for 4-site and 16-site 
clusters at 15% hole-doping are shown in Fig. [31 For 
both clusters and in the temperature range available, t' 
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FIG. 4: (color online) ^-dependence of T c at a) 5% and b) 15% 
dopings for different values of the dimensionless EP coupling, 
A (Eq. [2|. A finite t' increases T c in the under-doped region 
when EP coupling is present. For A = 0.6 in (a), the t' — 
point has been extrapolated form the data at temperatures 
larger than T c . 



suppresses the d-wave pairing susceptibility. 

In the presence of phonons, long-range hoppings change 
the polaronic effects which have a direct influence on T c . 
In previous works, we have found that at the intermedi- 
ate EP coupling, local phonons suppress T c due to po- 
laronic effects which reduce the mobility of carriers^. 
The polaronic effects are enhanced by the AF correla- 
tions 4 ^. Therefore, the effect of t' on the AF correlations 
will directly influence them. Using exact diagonalization 
methods, T. Tohyama el. a/. 10 ' 35 have shown that if sup- 
presses the AF correlations in the hole-doped cuprates. 
As a result, the polaronic effects are reduced by t' which 
enhances T c . 

When Holstein phonons are considered, T c increases 
with if in the under-doped region and remains almost 
unchanged in the intermediate-doped region as shown in 
Fig. 3J Note that without phonons, T c decreases with 
t' in the intermediate-doped and does not change in the 
under-doped region. The effect of if on T c becomes more 
significant with increasing A. As shown in Fig. 3] (a), for 
A = 0.6 at 5% doping, T c is strongly enhanced when \t'\ 
increases. At 15% doping, it is difficult to fix the filling- 
due to the large charge fluctuations for the values of A 
larger than 0.5. Presumably, T c would increase with t' 
at 15% doping for A = 0.6. Note that phonons do not 
increase T c , they only reverse the behavior of T c with t' . 
At fixed t', their effect is to reduce T c , but the reduction 
is less significant when |t'| is larger. 

Phonons change the behavior of P d Q with respect to 
if at low temperatures. As shown in Fig. [2] (b) and (d) 
with full symbols, t' has a small effect on suppressing 
Pdo around T c when phonons are present. While the 
band renormalization effects, caused by t', tend to sup- 
press Pdo , the reduction in the polaronic effects due to t' 
enhances P d Q. As a result of these two competing effects, 
Pdo remains almost unchanged near T c by changing t . 
On the other hand, the ^'-dependence of Vd is not influ- 
enced much by phonons. Therefore, T c increases in the 
under-doped region where V d shows a slight increase with 
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t'. 

It is known that the 2x2 cluster overestimates the d- 
wave superconductivity due to the neglect of phase fluc- 
tuations^. However, here we focus mainly on investigat- 
ing the relative dependence of T c on different parameters 
such as EP coupling and long-range hoppings and not 
on calculating the exact value of T c . The 16-site cluster 
results (shown in Fig. [3]) suggest that these trends do not 
change when larger clusters are considered. 

The transition temperature at small doping would be 
the strongest affected by phase fluctuations^. In gen- 
eral, t' should enhance T c by suppressing the phase fluc- 
tuations. In part, this is due to the suppression of AF 
correlations. The AF correlations reduce the mobility of 
the carriers and increase the effective mass which leads 
to the enhancement of phase fluctuations^,. This effect 
is enhanced in the presence of phonons. Phonons play a 
role similar to the AF correlations in enhancing the phase 
fluctuations by reducing the mobility of electrons. In ad- 
dition, through polaronic effects, they also enhance the 
AF correlations. Therefore, the effect of phonons on sup- 
pressing T c is underestimated in the absence of the phase 
fluctuations. The decrease of the polaronic effects due to 
t' would increase T c more significantly in the under-doped 
region in the presence of phase fluctuations. Hence, the 
effect of t 1 on T c is underestimated by small cluster cal- 



culations, suggesting that t 1 will have a stronger effect on 
T c in larger clusters. 

Conclusion- We find that without phonons, the long- 
range hoppings generally suppress T c in the Hubbard 
model. However, by including Holstcin phonons, T c in- 
creases with t' in the under-doped region for the hole- 
doped system while the suppression in T c due to t 1 is re- 
duced in the intermediate-doped region. Phonons do not 
increase T c , but rather reverse the behavior of T c with 
t'. We find that the increase in T c with t' becomes more 
significant for larger values of the EP coupling. We inter- 
pret these by the effect of t' on suppressing the polaronic 
effects as a result of suppressing the AF correlations and 
the interplay between AF and EP coupling. 
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